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A detailed description of Moscow-type (M-type) potential models for the NN interac- 
tion is given. The microscopic foundation of these models, which appear as a consequence 
of the composite quark structure of nucleons, is discussed. M-type models are shown to 
arise naturally in a coupled channel approach when compound or bag-like six-quark states, 
strongly coupled to the NN channel, are eliminated from the complete multiquark wave 
function. The role of the deep-lying bound states that appear in these models is elucidated. 
^ ' By introducing additional conditions of orthogonality to these compound six-quark states, 

a continuous series of almost on-shell equivalent nonlocal interaction models, characterized 
by a strong reduction or full absence of a local repulsive core (M-type models) , is generated. 
The predictions of these interaction models for 3N systems are analyzed in detail. It is shown 
that M-type models give, under certain conditions, a stronger binding of the 3N system than 
the original phase-equivalent model with nodeless wave functions. An analysis of the 3N 
system with the new versions of the Moscow NN potential describing also the higher even 
partial waves is presented. Large deviations from conventional NN force models are found 
' for the momentum distribution in the high momentum region. In particular, the Coulomb 

, displacement energy Eb( 3 He) — Eb( 3 H) - when Eb( 3 H) is extrapolated to the experimental 

value displays a promising agreement with experiment: AEc —740 KeV. The validity and 
limits of two-body NN potentials in nuclei is discussed in the light of our analysis. 
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I. INTRODUCTION 



It has been known for a long time that standard NN interaction models based on one-meson exchange 
jl]-[5|], meet significant difficulties in explaining standard properties of few-nucleon systems, ordinary nuclei, 
as well as nuclear matter In a number of cases, difficulties arise even in the well investigated area of 

low energies (~ 3 — 10 MeV) [^| [s]. In order to reach agreement between theory and experiment for some 
key observables (binding energy, magnetic moments, etc.) 37V-forces and meson-exchange currents (MEC) 
are very often introduced. However, in practical applications, the values of cut-off masses are sometimes 
incompatible with those used in the initial two-body interaction model Q . 

These discrepancies have resulted in a revival of interest in nonlocal nuclear force models (see e.g. PJlO|). 
The latter give predictions closer to experiment than those of conventional local models ||. In fact, from 
the viewpoint of the more fundamental level of quantum chromodynamics (QCD) and the quark model, 
the NN interaction must be strongly nonlocal at small distances |y],(l2| due to quark exchanges between 
the nucleons (three-quark clusters) |Il|dT3j|. A number of nonlocal models for the NN interaction, either 
phenomenological Jf4|,[l5| or microscopic ]f l|JT3] , ^6| 4l9[ , including both meson and quark degrees of freedom 
have already been suggested. However, most of them lead to rather complicated and bulky momentum- 
dependent NN interaction operators that can hardly be used for practical calculations in nuclear many-body 
systems. This explains the lack of qualitative studies of the structure of few-nucleon systems using quark 
degrees of freedom. Such studies would be extremely important for nuclear physics in general. 

Therefore, it is a rather urgent to develop a simple NN-interaction model which simultaneously takes 
into account both meson and quark degrees of freedom, and which is applicable to few- and many-nucleon 
systems. Having at our disposal a reasonably simple and realistic NN interaction model of the above hybrid 
type, we can study its qualitative and quantitative consequences in the many-body problem. 

The purpose of this work is to do just that. At first, we develop a simple nonlocal NN interaction model 
that is based on the quark structure of the nucleon. We then study its predictions for the three-body 
bound state problem and compare these predictions with those of conventional force models. An additional 
argument in favor of our approach is that quark degrees of freedom are not explicitly seen in the energy 
region {E™ <300 MeV). Only rather indirect signals of these hidden degrees of freedom can be studied. 
One such indirect quark effect considered in the present work is the internal fermization of the nuclear wave 
functions p0|. Such an "internal" fermization changes the wave function of the system mainly, but not 
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exclusively, at small distances and leads to the occurrence of additional internal nodes in the wave functions 
at about the same place where the repulsive core of conventional force models is situated. 

From a formal point of view, the elimination of quark degrees of freedom results in additional orthogonality 
conditions for the NN channel in conformity with the classical Feshbach method for projecting onto mutually 
orthogonal channels | 2l | . As a result, the orthogonality constraints lead to additional internal nodes in the 
wave functions p2| , [23| ]. Therefore, the height of the local repulsive core can be strongly reduced (see Sections 
III and IV) . We emphasize that a large part of the repulsion in lowest partial waves is not due to the local core 
produced by vector meson {uj and p) exchange but comes from these additional orthogonality conditions. As 
will become clearer later on, this leads to a substantial increase of the average nucleon kinetic and potential 
energy, which in turn has numerous consequences for the structure of the nucleus. Below we point out only 
the most important differences between conventional nuclear force models and M-type models. 

• Because of the additional orthogonality constraints, the strength of the repulsive core, mainly due to 
w-meson exchange, can be considerably reduced down to values, dictated by SU(3) symmetry i.e. 
SuNn/^ ~ 5. As a result, the attractive part of the potential is considerably deeper than that of 
conventional force models. However, this does not cause any discrepancy with experiment because the 
additional orthogonality conditions lead to an effective reduction of the local attraction at intermediate 
distances. 

• Due to the sharp increase of kinetic energies in the NN channel at small distances r^N < 1 fm, the 
relative importance of nonlocal and energy-dependent terms in the NN interaction operator coming 
from meson retardation effects Q or 6g-bags are strongly suppressed. As a result, all complicated 
and nonstatic short-range nonlocalities originating from different sources are replaced by a very simple 
nonlocal separable potential. This nonlocal separable potential has the form of a projection operator 
and effectively describes the additional orthogonality conditions [Q . This allows to drastically simplify 
the description of the NN interaction, especially at small internucleon distances. 

• An important difference is also the local character of the deep attractive potential which is universal 
for all partial waves. This universality can be interpreted as a consequence of pseudoscalar tt, rj and 
scalar tr-exchanges between nucleons. Here, the scalar tr-exchange, gives a strong attraction (depth 
~500 MeV p[ ) between nucleons at intermediate distances. 

For sake of brevity, we shall call NN- force models, which are characterized by (i) a strongly reduced local 
repulsive core (or even by the full absence of it), (ii) a strong attraction in the intermediate range and 
additional orthogonality conditions, as Moscow- (or M-) type models. For a substantiation of these models 
it is important to recall the following relatively new results. 

In a number of recent studies |^5|,g6| it was established that a limiting case of such models, which does 
not include any repulsive core at all, nearly coincides with the supersymmetric (SUSY) partner of the Reid- 
Soft-Core (RSC) potential. Thus, both types of interaction models, conventional force models with a local 
repulsive core and M-type force models, are related by purely algebraic symmetry transformations. For a 
pedagogical discussion of this point see Ref. Q . 

Recently Nakaishi-Maeda [p8| has found that the phenomenological separable NN-potential of Tabakin p^] , 
leading to internal nodes in the radial dcutcron and NN-scattering wave functions, is a simple unitary-pole 
approximation (UPA) to the above Moscow NN-potential, i.e. in the UPA-approach both models coincide 
with high accuracy. 

These results, established by independent groups, are not accidental; rather they are indicative of the 
quark dynamics underlying the Moscow NN-potential. One of the main purposes of the present work is to 
explain this in greater detail and to study the qualitative behavior of M-type NN interaction models in the 
nuclear three-body problem. 

This work is organized as follows. In Sect. 2 we provide a substantiation of M-type models within the 
framework of the quark model. Sect. 3 presents a pedagogicial introduction to the main physical ideas 
underlying the Moscow model. We study a simple scalar interaction and investigate its qualitative predictions 
for the three-nucleon bound state. Sect. 4 discusses the more complete model including tensor forces for 
positive parity partial waves. In Sect. 5 these more realistic models are applied to the three-nucleon bound- 
state problem. We emphasize the important contribution of the tensor force to the 3N binding energy. The 
interference of central and tensor forces inherent in these models is quite different from conventional force 
models. The differences in predictions between traditional force models and M-type models are pointed out. 
In the conclusion our main findings are summarized. 



II. ORIGIN OF NN INTERACTION MODELS WITH ADDITIONAL ORTHOGONALITY 

CONDITIONS: M-TYPE MODELS 
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A. Six-quark permutational symmetries and the NN interaction 



Consider a two-nucleon system that is described in the framework of the nonrelativistic quark model 
The quark Hamiltonian H consists of the following components: 

H = H + VoGE + V con f + VoME, (1) 

where Hq = 53j=i ( m » +Pi / (2wii)) is the kinetic energy of the 6g-system. The interaction Voge accounts for 
one-gluon exchange, and the confinement interaction V con f takes the standard form |]ll| , ^3| . In the one-meson 
exchange potential between quarks Vome we include, as usual, tt and u-meson exchange potentials [[l3 18 . 

The totally antisymmetric (denoted by subscript A) wave function ^a(x±, ...,xg) can be expanded into a 
sum of irreducible representations (IR) of the symmetry group Sq ■ The sum extends over all allowed 

IR of the symmetry group Sq, characterized by Young schemes [/], that lie in the outer product space S3 (3S3: 

y A (x u ...,x 6 ) = _L^*([/]r;ri, • r 6 )xcST ([/>), (2) 



nf lf]r 

where the coordinates x\ ■ ■ ■ x% collectively stand for the position, and the spin, isospin, and color quantum 
numbers of the six-quarks. Here, $([/]r; fi, r§) is the orbital part of the total six-quark wave function, and 
XcsT([.f]f) is its color-spin- isospin part. The Young scheme of the space part of the six-quark wave function 
is denoted by [f]r, where r is the corresponding Yamanouchi symbol. The Young scheme [f]f stands for the 
unique irreducible representation in spin-isospin-color space which is adjoint to [f]r. For brevity we omit 
the Young symbols [f)c, [f]s, [f]r, relating separately to the color, spin and isospin parts of the function 
Xcst- Furthermore, n/ is the dimension of the IR corresponding to the Young scheme [/]. 

One can easily see that, according to the Littlewood theorem [[H], the allowed Young schemes, which 
describe the permutational symmetry of six-quark orbital (X) wave functions $([/]r), are: 

[3]x x [3}x = [6] x + [42]* + [51] x + [33] x , (3) 

where the first two terms correspond to even orbital angular momenta (with positive parity), and the last 
two terms to odd orbital angular momenta (with negative parity) of the relative NN wave function. Thus, 
for even partial waves in the relative NN wave function, the allowed space symmetries are [6]x or [42] x, 
whereas for odd ones, they are \51]x or [33] x- Using the two-center shell model basis with distance R between 
the two centers, one can show 132] that all allowed six-quark wave functions of the type |S+S 3 [f]LST > 
approach the usual shell model configurations \s m p n [f]LST > (m + n = 6) in the limit R — ► 0. Thus, 
the NN state with the totally symmetric Young scheme [f]x = [6]x corresponds to the six-quark state 
|s 6 [6]x[2 3 ]c , s^ = 0,ST > while the NN state with mixed symmetry configuration [42] x corresponds to a 
six-quark state with two excited p-shell quarks: 

\s i p 2 [42]x[f}csL = 0,ST>, 

where [f]cs = [42], [321], [2 3 ], [31 3 ], [24 4 ] are the possible Young schemes appearing in the Clebsch-Gordan 
decomposition of the internal product [2 3 ]c[42]s. 

In a series of works by many different authors ]^3|^6|j33| -pq] , it was shown that, with a quark Hamiltonian 
of the form of cq.(l), the excited six-quark configurations |(0s) 5 (ls)[42] x > and |(0s) 4 (0p) 2 [42] x > which 
are compatible with S-wave relative NN motion, are admixed with large weights to the fully symmetric 
six-quark configuration |(0s) 6 [6]js s : >. The color-magnetic term ~ • \j<Ji ■ <Tj is mainly responsible for 
this. Thus, in lowest partial waves of the relative NN wave function, there is a superposition of two different 
6g-space symmetries: the fully symmetric | (Os) 6 [6] > and the mixed symmetric |(0s) 4 (0p) 2 [42]x > for 
S'-waves, and similarly the |(0s) 5 (0p)[51]x > and |(0s) 3 (0p) 3 [33]x > for P-waves. It is important to remark 
that the above superpositions include the excited p3 single-quark states. 

In recent work |l3| ] , it was shown that six-quark components with different space symmetries play a very 
different role in the NN interaction. For example, the totally symmetric six-quark components |(0s) 6 [6]x > 
are projected onto the NN cluster channel with rather small weights, whereas they have large projections 
onto the AA and hidden-colour channels. In contrast, the mixed symmetry components |s 4 p 2 [42]jf > (in 
S'-waves) have large projections onto both the cluster NN channel and the nucleon-isobar channels N^N^- 

Thus, there is a natural separation of the complete six-quark wave function into two mutually orthogonal 
parts of a different physical nature: 

®A = *6<z + Vclust, (4) 

where ^e q is the bag-like component, which can be constructed from square integrable functions tpi corre- 
sponding to various six-quark bag states. These are the states with maximal space symmetry, i.e. s 6 for 
relative S'-waves and s 5 p for relative P-waves: 
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^eg^^C^, (5) 

i=l 

where N is the number of six-quark states. For the cluster component we use the standard RGM-ansatz: 

*d««t(ei,e 3 ,R) = ^{^(CiVjv(e 2 )x(R)}, (6) 

where A is the antisymmetrizer with respect to the six quarks. Here, <v9jv(£j) is the three-quark wave function 
of a single nucleon where the internal coordinates are collectively denoted by ^; and x(R) is the relative 
motion wave function of two three-quark clusters. 

B. Two-phase model for the NN interaction 

Our approach differs from the majority of hybrid NN interaction models which are also based on a 
decomposition of the form of eq.(4). In contrast to, for example, the quark compound bag (QCB) model by 
Simonov [14], we require the mutual orthogonality of the components and ^ c i us t'- 

< ^6q\^clust >= 0. (7) 

Moreover, we require that the cluster component be orthogonal to all bag-like states ifi from which the 
component ^>Q q is constructed: 

< ViWclust >= 0, i=l,...N. (8) 

When combined with eqs. (4)-(6), this leads to corresponding orthogonality conditions for the RGM relative 
motion function x{R): 

< 9t \x>=0, i=l,...N, (9) 

where 

ff< (R) =< Pifo.ea.BOIMfiWk) > ( 10 ) 

are projections of six-quark bag-like functions onto the TV TV-channel. To emphasize the importance of 
orthogonality conditions in our approach we put a ~ sign over the relative motion wave function \- We 
point out that a similar model with mutually orthogonal 6q and NN channels in configuration space has been 
suggested by Lomon |jl9|| . The model of Lomon is similar in spirit to our model but different in realization. 

It is obvious that at intermediate and large distances, where the cluster-like components dominate, the 
NN dynamics should be rather well described by a meson-exchange (one-boson-exchange(OBE) and two- 
pion exchange (TPE)) model. Explicit quark and gluon degrees of freedom are unimportant here. On 
the contrary, in the short-distance regime, with maximal overlap of the nucleon wave functions, explicit 
quark-gluon degrees of freedom, described by the compound states <pi should play a decisive role. In this 
way we arrive at what could be called a duality principle for the baryon-baryon interaction which can be 
formulated as follows. The total six-quark wave function naturally separates into two mutually orthogonal 
(i.e. nonoverlapping) components: 

• One component has a three-quark cluster structure where the baryons (including isobars) can be 
considered as separate entities, and where the dominant dynamical mechanism is meson exchange 
between the clusters, while the quark-gluon degrees of freedom are only of minor significance. 

• The other component has a six-quark bag structure. Its dynamics is governed by explicit quark-gluon 
degrees of freedom. This component is sensitive to for example, the particular form of confinement, the 
value of the scalar quark condensate < qq >, etc., and depends only weakly on the "external" meson 
dynamics. 

Therefore, when constructing a proper NN interaction potential, the bag-like components should be some- 
how excluded from the very beginning, because these components are very hard to describe by any reason- 
able NN potential. For this component one should use a different formalism. This is in complete analogy 
to the optical potential model in nuclear physics, which is not at all applicable in situations, where the 
nucleon-nucleus scattering proceeds via isolated compound-states of the nucleus. In this analogy the op- 
tical nucleon-nucleus potential corresponds to the meson exchange dynamics between nucleons while the 
nucleon-nucleus compound states correspond to the six-quark bags in NN scattering. 
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Projecting the total six-quark Schrodinger equation 

= 

onto the six-quark states ifi and onto the NN cluster channel we obtain, using the orthogonality conditions 
(8), the following set of coupled equations: 

/ £.•(#« - es 13 )c 3 + < Mx > = o, (ii) 1 



£j Qfj + (fl RGM - EJ\f)x = 0- (12) 



Here, we have introduced the following abbreviations 

=<<Pi\H\<pi>, (13) 

/i(R) =< <fi\H\<p N <p N > . (14) 
-j^rgm an( j _yyjr are integral operators with kernels: 

H RGM {Il,IL')=<<p N v N \HA\vNVN>, (15) 

j\f(R, R') =< (pN ( fiN\A\(fiNfN > ■ (16) 
Solving the linear algebraic equations (11) with respect to the coefficients Cj and substituting its solution 

C i = Y i (E-H)^<f i \x> (17) 

3 

into eq. (12), we find the equation of motion for the NN relative wave function X'- 

U RGM - EAf+J^lfi > (E ~ H)^ 1 <fi\\x = (18) 
v y / 

This equation must be solved together with the additional orthogonality conditions (9). After that one can 
calculate the coefficients Cj via eq. (17) and hence the total wave function Eq. (18) constitutes our 
two-phase model for the NN interaction. 

Due to the appearance of the norm kernel Af, the RGM relative motion wave function x i n ec l- (18) 
cannot yet be interpreted as the probability amplitude for finding the nucleons at a relative distance R. In 
order to obtain a NN relative motion wave function that can be interpreted in the usual sense, we define a 
renormalized function x — ■N' 1 ^ 2 X an d a corresponding effective nucleon-nucleon interaction: 

(T fl + V eff - E) x = 0, (19) 

where 

V eff = jV'^H^jV 1 ' 2 -T r + J2 \fi >(E- H)^ 1 < fa, (20) 



The effect of the overlap kernel 7V(R, R') for the NN system can be reasonably well approximated by the 
expression 

AA(R,R')^<S(R-R), (21) 

which differs from the usual local NN potential just by a combinatorial factor. Thus, the effective interaction 
can be cast into the form: 

^e// = Vdir + V ex + Vff q N , (22) 

where Vdir is the direct (folding) potential, including one-meson exchanges between three-quark clusters 
only, i.e. the one-meson-exchange nucleon-nucleon potential. Because the one-gluon exchange operator is 
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diagonal with respect to quark permutations between nucleons it does not contribute to Vdir- V ex is the 
nonlocal short-range exchange potential, and VNqN is an effective potential due to the coupling of NN- and 
6c/-channels: 

V NqN = 10 £ \h >(E- H)^ 1 < f s \. (23) 

ij 

The effective two-nucleon Schroderinger equation for the orthogonalized and renormalized relative motion 
wave function x(R): 

{T R + V dir + V ex + V NqN -E)jc = , , 

< 5l | X > = 0, i = 1,...N { ' 

provides the basis for the justification of NN-interaction models of Moscow type. 

The main point is that the solution of eq.(24) is defined only in the subspace orthogonal to the functions <?j. 
Due to this orthogonality requirement, we inevitably introduce nodes in the NN scattering wave functions at 
small distances [^2|,^3],[33] |3(| . The node positions are very stable as the NN scattering energy E is increased 
jjJLpSi. This behavior of scattering wave functions can be quite naturally formulated in the language of 
interaction potentials that have a very deep attractive well (~1 GeV) pj| with additional deep- lying bound 
states. In fact, compared to the deep attractive well with \Vq\ S3 1 GeV, the relative scattering energy E is 
negligible. Therefore, the position of the innermost node r no de is basically determined by the well depth Vo: 

y/m N (E — V Q ) r node S3 ^Jm N \ V \r node S3 ir. 

Thus, the inner node of the scattering wave function is almost stationary for center of mass energies up to 
~2 GeV in the laboratory frame in nonrelativistic kinematics. Inserting the numerical values for the nucleon 
mass, and Vo we find that the stationary node is located at R n ~0.6 fm in the case of relative S'-waves 
[R n ~0.9 fm in the case of P- waves). We point out that this node plays the role of the repulsive core in the 
NN interaction. 

Because the nodal solutions of the basic equation (24) with orthogonality conditions correspond to very 
large kinetic energies, the short-range nonlocal potential V ex , which includes the quark exchange diagrams, 
plays only a minor role in the subspace of nodal wave functions. This is in analogy to the nuclear Orthog- 
onality Condition Model (OCM) of Saito |3(|. Therefore, V ex can be omitted from equation (24) in good 
approximation. 

Finally, we take the last step in the substantiation of the model. In absence of the vector (p and lu) meson 
contribution in the total Hamiltonian (1), the meson-exchange potential Vdir in eq. (24) corresponds to tt 
and cr-meson exchange. Therefore, it is strongly attractive with a depth of the central potential of about 
^500 McV. Together with the nonlocal attractive (for E < E\) part VnoN , it results in an effective deep 
(~1 GeV) attractive potential, which was found in our initial attempts p3|, and has since been justified 



||13|,[22 J27| 37j] within the framework of the quark model. In fact, we have found out only quite recently []38| 
that the combination of the local attractive cr-meson exchange potential (with a depth of about 500 MeV) 
and the nonlocal attraction VNqN, representing the coupling of the NN and 6q-bag-like channels, is exactly 
phase-shift equivalent to the deep (~1 GeV) effective attractive potential of M-type. 

The Moscow potential is the simplest local model, which ensures the orthogonality of its solutions to a 
function ipo that approximately describes the projection of the fully symmetric (0s) 6 six-quark state onto the 
NN-channcl (in S'-waves) . We emphasize that the occurrence of unobscrvablc deep- lying bound states in the 
equation (24) poses no problem, because the equation for the NN channel wave function is considered in the 
appropriate subspace of nodal wave functions, which is orthogonal to the nodeless bound states. Additional 
extra deeply-bound states arise only in the solution of the appropriate Schrodinger equation in the complete 
space. 

In order to understand the reasons for the occurrence of deep-lying bound states, it is very instructive 
to recollect [ fL3| that, the weight of the NN-configuration (with unexcited nucleons) in symmetric bag- like 
six-quark states | s 6 [6] [2 3 ]r;sl/ = 0, ST > is very small; whereas the weight of the AA channel is appreciably 
larger and the main contribution comes from hidden colour CC channels, whose wave functions are localized 
(because of confinement) at small distances |3!|. As a result of the strong coupling of the initial NN channel 
and the bag-like six-quark channels, the probability for finding an NN configuration with unexcited nucleons 
vanishes almost completely at small distances^]. Thus, at very small distances, where the extra deeply-bound 



1 A similar conclusion concerning the suppression of the NN-channel in favor of the AA- and other isobar channels 
can possibly be derived in a formalism based on meson Lagrangians |40| without any reference to the quark model. 
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states are localized, there is basically no nucleon-nucleon configuration. In this regime, the 2iV-system 
can be viewed as though "being dissolved" in a quark "soup". Remembering the Cheshire Cat smile, it is 
appropriate to name these deep-lying bound states: "Cheshire-Cat Bound States (CCBS)". 

In applications of the model to many-nucleon systems, it should be kept in mind that the Moscow potential 
describes in an effective way the single- channel system, i.e. the 2N system with unexcited nucleons, while 
the real 2N-system should be described by the above two-component model. Thus the description of three- 
nucleon systems given by the Moscow potential model must be supplemented with an additional 6g+N 
contribution which is expected to be strongly attractive. Before we discuss the Moscow NN potential in 
greater detail, we explain the main ideas of our approach within the framework of a simple toy model. 

III. A SIMPLE M-TYPE MODEL FOR THE NN- INTERACTION 

A. General considerations 

In this section we compare two alternative models for the NN-interaction. For clarity, we use the simplest 
possible example, i.e. two spinless nucleons interacting via a scalar potential V(r): 

V(r) = g R V R (r) + g A V A (r). (25) 

The potential includes a short-range repulsive core gnVji(r), and a long-range attraction g A V A (r). The 
coupling constants g R and g A are chosen such that the system has a single bound state with energy Eq 

("deuteron") and scattering phase shifts Si(E), Z=0,1, 

Let us try to find a modified potential V'(r) — g' nVji(r) + g A V A (r) with a reduced strength of the 
repulsive core, but simultaneously introduce an additional constraint that ensures the orthogonality of the 
solution of the Schrodinger equation to some function ip localized in the core area: 

f (T r + V'(r) - e)f(r) = (26a) 
\ (<Po\fc)) = (266) 

The physical meaning of the additional orthogonality constraint has already been discussed and partly been 
elucidated above. Bere, we consider this orthogonality condition as a formal constraint that the solution has 
to satisfy. 

Equation (26a) with the additional condition (26b) can be rewritten in the form of a single equation with 
the nonlocal interaction (see [j36f): 

(T r + V'(r) - e)$(r)) = ko><¥>o|T r + V'(r)$(r)). (27) 

The solutions of this equation are, as one can easily see, orthogonal to ipo for any e =/= 0. Only for e=0, 
the required orthogonality is not guaranteed and Eq.(27) gives an incorrect behavior for the corresponding 
off-shell t-matrix. 

One can apply a more convenient approach, known as the method of orthogonalizing pseudo-potentials 
(OPP) J24|. In this approach, the additional orthogonality condition (26b) is taken into account via the 
projection operator T = \<fo)(<fo\ which projects onto the "forbidden" subspace with a large coupling constant 

(T r + V'(r)+ fJ ,\<p Q )(<p Q \-e)\i> fi (r))^ oo = 0. (28) 

In our example, the forbidden subspace is just a one-dimensional subspace spanned by the vector ipQ. In the 
limit /i — * oo, the solution of equation (28) ^(r) have been shown [|24f to be rigorously orthogonal to ipo, 
i.e.: 

lim (ip \i> l _ l ) = 0. 

fi — >oo 

For £ / the solution i/j^ coincides with the solution of the Saito equation (27). The projection operator 
r = l^oX^ol defines a "forbidden" subspace Hr of the full two-particle Bilbert space Ti which can be 
decomposed as 

7~ir (J) Hq = H 

where H.q is the orthogonal complement to 7ir- 

Thus, we consider the problem of finding a modified potential V'(r) that is phase-shift equivalent to the 
initial potential V(r) of eq.(25) and which has the bound state at the same energy e in the subspace Hq 
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orthogonal to ipo ■ For an arbitrary ipa , rigorous solutions of such an inverse scattering problem are, to our 
knowledge, not available. Therefore, we search for approximate solutions. In order to further elucidate the 
relation of M-type models and conventional NN force models, we are interested in such solutions that contain 
a repulsive core that is weaker than that in the initial potential V(r). Although the modified potential V'(r) 
becomes deeper when the core is weakened, the introduction of the additional orthogonality constraint of 
eq. (26b) renders the Hamiltonian effectively weaker. Due to this compensation, both the phase shifts and 
the bound state energies remain the same, albeit with some error. 

On the other hand, if we use the OPP-method |24|,|3(| the modified pseudo-potential has the form: 

V'(r) = V'(r) + fx\cp )(<p \ (29) 

Using this method, it is possible to work in the total space instead of only in the subspace. This is one of 
the practical advantages of the OPP-method. Thus, in the OPP approach, a part of the local short-range 
repulsion is replaced by a separable repulsive potential with an infinite coupling constant /i (29). If (po is a 
bound-state solution of the potential V'(r), the potentials V' and V' are completely phase-shift equivalent. 
For this case, the solution of the inverse problem, i.e. the transition to a phase-shift equivalent potential 
with a local repulsive core without a deep-lying bound state, is well known. This is a supersymmctrical 
(SUSY) transformation [^5|,|6). However, when extra bound states (number n) are eliminated, a SUSY- 
transformation always introduces a repulsive core of centrifugal type, i.e. v s (r) ~ n(n + l)/r 2 , instead 

of the usual Yukawa-type core in local OBEP models. Furthermore, the invariance of phase shifts under a 
SUSY transformation holds only for a given partial wave. 

Here, we want to consider the more general case of a repulsive core of intermediate strength. Furthermore, 
we want to preserve the invariance of the phase shifts in several not just in one partial waves. Most 
importantly, we do not want to introduce any superfluous bound states into our model from the very 
beginning. Therefore, we refrain from using the rigorous results of SUSY-transformations and we look for 
a potential that is close to the phase-equivalent potential in the subspace H.q orthogonal to a given bound 
state. 



B. The Malfliet-Tjon potential 

To be specific let us consider a particular example. As the initial potential we take the Malfliet-Tjon 
potential MT-V |4l) . The MT-V potential is a simple scalar model for the NN interaction, which yields an 
average deuteron binding energy of —0.4136 MeV and describes the average 3 Si and 1 S'o NN-phase shifts 
reasonably well up to energies of about 300 MeV: 

t t , s exp(-^ fl r) cxp(-^ A r) 

V M T-v(r)=g R h g A , (30) 

where the parameters values are 

g R = 1458.05 MeV, g A = -578.09 MeV, \i R = 3.11 fm, [i A = 1.55 fm. (30a) 
Next, we introduce an orthogonality condition to the function <po( r ) of Gaussian form: 

^ (r) = Vexp(-^. (31) 

This (0s) harmonic oscillator Jh.oA function approximates quite well the projection of the six-quark bag 
| s 6 [6]) onto the NN channel Jl3| , p3| , p7][39[ . For various values of the h.o. radius ro we now determine the 
coupling constants g' R and g' A of the modified potential: 

, , cxp(-^ fl r) , exp(-^r) 

V r = 9r + 9a (32) 

(J>Rr fJ. A r 

acting in the subspace Ti.Q, which ensure the best fit to the S- wave phase shift for the initial potential MT-V 
(30), acting in the full space Ti. For the sake of simplicity we do not change the inverse radii of the core, 
Hr, and of the attractive potential \x A . 

The results of this fit for different values of radius ro are presented in Table 1. It can be seen that a fit is 
possible for ro lying in the interval 0.15-^0.5 fm. Thus, in the subspace TLq orthogonal to ¥?o( r )> the partial 
,5- wave phase shifts for the potential V'(r) of eq.(32) reproduce fairly well the corresponding phase shifts of 
the initial potential (30) in the energy range -j- 300 MeV (see the % 2 values in Table 1). We do not present 
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here the respective figures because the difference between the phase shifts of the initial potential MT-V and 
those of V'(r) is almost indistinguishable in the graph. 

It should be noted that the local core disappears already at fm. A further increase in ro results 

in an increase of the radius r no d e of the internal node in the wave function and of the amplitude of internal 
loop specified by the ratio (3 of absolute values of the wave function in the first and in second maxima (see 
the fifth column in the Table 1). After such a modification of the potential, the bound state energy £o (in 
the full space) varies from —2.51 up to —600 MeV. However, if we require the orthogonality to ipo, i.e. if we 
are operating in the orthogonal subspace, the bound-state energy e"o in Table 1 remains almost constant that 
is of the order of 1 MeV. 

Thus, we obtain a series of approximately spectral-equivalent Hamiltonians^j. The best (i.e. the small- 
est)value of x 2 is obtained for ro=0.4 fm when there are two bound states in the potential V'(r): the first 
one is deep-lying with the energy eo = —528 MeV (see the third row in Table 1) and the second one is 
close to the NN-threshold at the energy e\ = —0.4 MeV, the eigenfunction of the deep bound state being 
very close to the h.o. function fo{ r )- in other words, we arrive again at a deep attractive potential without 
any local repulsive core but with an eigenstate projector, that is, to a construction that is very similar to 
the Moscow potential model. Thus, there are two extreme cases: (i) the deep purely attractive potential 
with the "extra" deep-lying state at eo = —528 MeV corresponding to ro=0.4 fm, and (ii) the initial MT-V 
potential which does not contain the "extra" deep- lying bound state state at eo = —528 MeV but has a local 
repulsive core instead. 

These two phase shift equivalent potentials are approximate supersymmctric partners. All other cases 
in Table 1 can be considered as intermediate variants, lying between these extremes^]. By changing the 
"projector radius" ro from zero to 0.4 fm, we perform a continuous transformation from one extreme model 
to the other. In doing so, the position r noc i e of internal node of wave function changes and the amplitude 
of the internal loop of wave function increases. This can be seen in Fig. 1, which shows the scattering wave 
functions at 1 MeV for three potentials: the initial local MT-V (dashed line 1), the intermediate nonlocal 
potential with ro=0.2 fm (solid curve 2) and the potential having an approximate eigen state projector with 
ro=0.4 fm (dotted curve 3). 

We emphasize again that, the second extreme case with ro=0.4 fm involves an almost exact eigen state 
projector, i.e. for the given values of the coupling constants g' R and g' A (see Table 1) the modified potential 
(32) leads to a deep-lying bound state with with energy eo — —528 MeV. The corresponding eigen function 
is very close to <fo(r) at ro=0.4 fm. Because the Hamiltonian is Hermitian, all scattering wave functions, 
and also the eigenfunction of the second (near-threshold) bound state at £\ = —0.41 MeV are automati- 
cally orthogonal to the wave function of the deep-lying level. Hence, there is no need for any additional 
orthogonality conditions. 

As a result, the supersymmetric partners are local potentials, while all phase-equivalent intermediate cases 
correspond to nonlocal interactions in the full space, or alternatively, to local interactions in the subspace 
TLq- However, according to our intention, we would like to obtain a good description not only in partial 
S'-waves, but also in D-waves and possibly in other even partial waves. In a realistic hybrid model of 
the baryon - baryon interaction, a quark bag is formed only in the lowest partial waves. Therefore, the 
additional orthogonality condition must be included only in lowest partial waves (S- and P-), whereas the 
same Hamiltonian without any additional orthogonality constraints is expected to correctly describe the 
phase shifts in higher partial waves. 

Thus, we attempt to find a modified potential (32) acting in the subspace TLq that simultaneously describes 
the S- and D-wave phase shifts for the initial MT-V potential. This is possible with acceptable accuracy for 
various values of 7'o. Table 2 presents several variants of the modified potential (32) with different values ro 
that are fitted to the S- and D-wave phase shifts under the condition that these are orthogonal to the S— 
wave function tpQ. Because we intend to use the above potentials in three-body calculations, we present here 
the versions for which the "deuteron" binding energies are practically identical (=0.413 MeV). The results 
of three-body calculations with these potentials is also given in the Table 2 and are discussed below. 

Figs. 2 and 3 show a comparison of the S- and D-wave phase shifts in the energy range 0-400 MeV for 
three potentials: (i) the initial local MT-V model, (ii) a modified nonlocal potential V'(r) with ro=0.2475 fm 
(i.e. a potential without a repulsive core that is close to the SUSY-partner of the MT-V), and (iii) a modified 
nonlocal potential V'(r) with ro=0.15 fm (an intermediate variant). It is clearly seen that the potential with 



2 In contrast to trivially spectral-equivalent Hamiltonians that are obtained as a result of the unitary transformation 
H' — U~ X HU and which have been investigated in detail in the literature [Q, our transformed potentials do not 
include any velocity dependence. Most importantly however, they have a different physical meaning. 

3 Such an intermediate model may be relevant for modelling the physical w-meson exchange interaction with a 
reduced wNN coupling strength. 
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the strongly reduced repulsive core acting in the subspace TLq provides a satisfactory description of both S- 
and D-wave phase shifts of the initial MT-V potential acting in the full space 7i. 

In this way, our simple toy-model shows that the introduction of an additional condition which ensures 
orthogonality to some localized state ipa enables us to perform a continuous transition between the local 
core model and the model with "extra" bound states but no core. Both alternatives describe the interaction 
between composite nucleons equally well. The extreme models happen to be SUSY partners. In between, 
one has a series of almost phase shift equivalent potentials, differing by the core strength and by the spatial 
extension (ro) of the state ifQ. 

C. Consequences for three-nucleon bound states 

Next we discuss, how the properties of the 3N-system change when going from one NN-force model to 
the other. In table 2 we compare the 3N predictions of the initial MT-V model (30, 30a) with those of the 
different variants of the orthogonalized model. Table 2 displays bound-state energies -E3, kinetic energies 
T3, and also the minimum and secondary maximum position of the 3 H charge form factor together with its 
value F max in the secondary maximum for different variants of the potential (32) and different functions (i.e. 
different values of ro) ¥>o( r )- 

One can see that an increase in the "projector radius" ro from up to 0.25 fm leads to an increase of the 
three-body binding energy from 8.25 up to 8.99 MeV. As ro increases further up to 0.4 fm the binding energy 
begins to decrease again. These results contradict the results of Nakaishi-Maeda who showed, that the 
replacement of the local core by the orthogonality condition leads always to a reduction of the three-particle 
bound-state energy. According to jl3|], the "stronger" the imposed orthogonality condition, i.e. the larger 
the amplitude of the inner loop in the radial wave functions, the more pronounced the reduction of the 
binding energy. At first sight, our results seem to contradict the 3N-calculations with the Moscow potential, 
undertaken by several groups (see e.g. |d|-|4(|). These calculations showed that the replacement of the 
local repulsive core NN-potential (for example, the Reid potential) by the Moscow AW-potential leads to a 
reduction of the binding energy from 7 down to 4.5 MeV fjij . 

The difference between our results and the Nakaishi-Maeda findings is explained as follows. The transition 
from the local core to the orthogonality condition causes a sharp increase in the kinetic energy of the 3N 
system (see the fifth column in Table 2). This in turn leads to an increase of the contributions of higher 
partial waves. The main role is played by an attraction in D-waves - see the seventh column in Table. 2. It 
should be emphasized that the MT-V potential is central and therefore does not mix three-body states with 
different values of the total angular momentum L. We consider here the mixing within a state of given total 
angular momentum L (e.g. L~ 0). By higher partial waves, we mean states with A = I = 2,4, where A 
and I are angular momenta associated with the Jacobi coordinates r and p, respectively. 

The change in the 3N binding energy thus depends on the magnitude and sign of the D-wave interaction: 
if this interaction is strongly attractive, the binding energy is increased, if it is repulsive the binding energy 
is decreased. To verify this statement, we have carried out 3N calculations with strong repulsion in .D-wave 
interactions, while keeping the S-w&ve interaction fixed. The results are presented in the eighth column 
of Table 2 {E^ ep ). With increasing ro, i.e. with a strengthening of the orthogonality constraints, the 
contribution of the internal D-wave grows and the binding energy E^ ep with repulsive D-wave interaction 
monotonically decreases. 

The calculations of Nakaishi-Maeda |l3| and Hahn et al. [|44| were based on the Faddeev equations and 
ignored usually small contribution of higher partial waves. Under these restrictions, our calculation also 
shows that the orthogonality condition always leads to a stronger effective repulsion than the original local 
repulsive core model.n As a result, without the higher partial waves the binding of the 3N system becomes 
weaker as one goes from the local repulsive core model to an M-type model. If we include the contribution 
of higher partial waves in pair subsystems (for a fixed value of the total orbital angular momentum L), the 
3N binding energy changes in the way described above. An appreciable increase of the contribution of the 
higher partial waves in pair subsystems due to the increase of the internal kinetic energy is the distinctive 
feature of M-type models. A more detailed discussion of this property is given below when we consider a 
more realistic model. 



4 This is probably due to the fact that the radius of the separable repulsive core in the pseudo-potential fi\ipo >< ipo\ 
(for n — » 00), is much greater than the radius of the original local repulsive core. 
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IV. THE MOSCOW POTENTIAL MODEL 



A. One-channel model 

Early attempts to construct an NN-potential with an additional "forbidden" state, i.e. a model of M- 
type, were undertaken already in the middle of the seventies |23f| . However, in these early works only 
central potentials and only even partial waves were considered. We point out that these early attempts were 
undertaken yet before the color degree of freedom quark was fully established. But for colourless quarks 
(fermions), the lowest (as it then seemed) six-quark configuration |(0s) 6 [6]js:L = 0, ST > is strictly forbidden 
by the Pauli principle. Therefore, these NN potentials were interpreted in complete analogy with the effective 
potentials between nuclear clusters (e.g. the a — a potential), in other words as potentials with bound states 
forbidden by the Pauli principle. 

The first version of a realistic Moscow potential was suggested in 1983 [^J and further details can be 
found in Ref. p8[ . This version describes the NN interaction only in the 1 Sq and 3 Si — 3 D± channels. Later 
on the description was extended to the lowest odd parity partial waves ( 3 Po, 3 Pi, 3 P2 — 3 P2) [^|- A new 
and improved version of the Moscow potential (version 86) was published in ]5C| . This version still had a 
node in the deuteron D-w&ve for which there is no evidence in microscopic quark models |u] , |l3| , [l6| . Finally, 
in 1990, we have changed the truncation of the tensor OPE potential p7| , which eliminated the D-w&ve 
node (version 90). 

The main difference between the Moscow potential and standard models for the NN interaction is the 
absence of a local repulsive core at small distances. Instead of the core, the interaction is described by a 
deep attractive potential Vq exp(— rjr 2 ) plus an appropriate condition of orthogonality to an extra deeply 
bound state. Owing to the deep-lying bound state, the physical NN wave functions have a node at small 
distances. The position of this node at around ~ 0.6 fm is almost stationary with increasing energy. In 
our approach, this stationary node plays the role of the repulsive core in standard models, and provides 
the correct behavior of the NN scattering phase shifts. The physical meaning of the "extra" bound states 
in the lowest partial waves has been repeatedly discussed in the literature ^3,p7],[37| . They either simulate 
six-quark bag compound states, for example, |s 6 [6]x) in S-waves, and |s 5 p[51jx)in relative P-waves, which 
cannot be adequately described in terms of nuclcon degrees of freedom only Jl3| , ^5| , pT|j2^ ] ; or they describe 
bound states in the nucleon plus Roper resonance NN*(1440) channel p7|| . We recall that the interaction 
in the NN*(1440) channel must be very similar to the original NN interaction because the Roper resonance 
has the same quantum numbers as the nucleon. 

Because we are mainly interested in the predictions of Moscow-type NN-potentials for 3N systems, a more 
complete version of such potentials including even partial waves with I <4 is necessary. Here, we present a 
few new versions for even partial waves. The form of the interaction for odd partial waves (for the version 
86) is discussed in J5(J. An updated version of the Moscow potential for odd partial waves will be considered 
in a future publication. In accordance with general expectations, six-quark bags appear mainly in the lowest 
S and P partial waves (for < 800 MeV). Therefore, the extra bound states describing the projections 

of these six-quark states onto the NN channel can occur only in these lowest partial waves. 

In the higher partial waves there are no extra bound states and correspondingly no nodes in the wave 
function. However, an interaction that has in the lowest partial wave only an attractive local potential and 
a repulsive projector multiplied by a large positive constant, is too strong for higher partial waves. In order 
to take into account the short-range repulsion generated by w-meson exchange, and in order to retain the 
universal form of the interaction in all partial waves, we add a repulsive core in the same separable form 
with a Gaussian form factor, as in the lowest partial waves, but with a finite positive coupling constant. We 
will then arrive at some intermediate model of the type discussed in the previous section. 

The final Moscow potential consists of two parts. First, a local ^-independent part, which explicitly includes 
the one-pion-exchange (OPE) potential and a deep attractive well Voexp(— rjr 2 ) that depends on the spin 
and parity of the NN system. Second, a state-dependent separable repulsive core with Gaussian form factor, 
which provides the correct P/-wave phase shifts. This form of the potential is not only universal, but it is 
also convenient for three-body variational calculations, because it is possible to calculate all of the matrix 
elements of the Hamiltonian analytically |[l] [33| . The separable core providing the short-range repulsion 
depends on I and J. Therefore, it is not necessary to explicitly include a spin-orbit interaction for even 
partial waves. 

The OPE potential is truncated in a suitable way at small distances (see below). The tensor interaction 
which couples partial waves with angular momenta I and I ± 2 is in all partial waves quite reasonably 
described by the truncated OPE-potential. Thus, the channel coupling is easily controlled by the truncation 
parameter. For partial waves with I > 4, the repulsive core is not required (for Ei a \, < 400 MeV) since the 
phase shifts with I >4 are completely determined by the long-range OPE tail of the interaction. Thus, the 
NN potential for even partial waves has the form: 

VnN = Vloc + V sep , (33) 



11 



where Vi oc is a local I- independent part of the interaction which, however, depends on spin and parity: 

V loc (r) = V c (r) + V T (r)S 12 ; (34) 

V c (r) = VoeM-Vr 2 ) + V° PE (r) ■ f tr (r); (35) 
V T (r) = V° PE (r)-(f tr (rT; (36) 



where the cut-off factor is 



ftr{r) = 1 - exp(-ar). (37) 



The exponent n is different for various versions of the model. The OPE potentials Vg rE and V^ ybl have 
the standard form: 

Vg PE {r) = V OPE exp{-x)/x; x = ^r; 

Vf PE (r) = V OPE {1 + 3/x + 3/x 2 ) exp(-x)/x. 
We use an averaged 7r-meson mass fj, and a it NN coupling constant g 2 NN /^ = 13.8 [ f^ , pl|]5"c| ] 

v ope = _?InnJI^ = „i . 69Me v ; M = 0.6995 fm- 1 . 

The separable repulsive core V sep has the form (for 1=2 and 3): 

V sep =\\<p><tp\; (38) 



tp = Nr l+1 exp {-) j ; J ip' 2 dr = 1. (39) 

The strength constants A for the separable core vary for different J and I, whereas the radius of the 
repulsive core tq varies only slightly from state to state (see Table 3). For I > 4 the term V sep is absent. 
Parameters for several variants of the iV TV-potential differing in the truncation of the tensor part in triplet- 
even channels are given in Table 3. We have found (see Sect. 4) that the properties of NN- and 3N-systems 
calculated with M-type potentials (this is also true for conventional NN potentials [|||rj ||,|54|] strongly depend 
on the behavior of the tensor potential at small distances. In the 1986 version of Moscow potential (variant 
(A) in Table 3), an exponential truncation of the tensor part (see eq.(37)) with a value n=3 and an inverse 
radius of truncation a = 3 fm -1 was used. In that case, the tensor potential truncated according to eqs. 
(36) and (37) corresponds to a large negative constant at the origin (Vt(0) = —5500 MeV). Therefore, there 
is a very strong coupling of the 3 Si — 3 Di channels at small distances. This coupling results in a large loop 
in the deuteron D-wave, as well as in a large D-state admixture Pd of 6.4%. The relative amplitude of this 
loop A E oop compared to the maximum value of the -D-wave function is given in Table 3. 

However, in view of the relatively weak quark-quark tensor force, it is desirable to have an interaction 
variant, for which the node in the deuteron D-wave is absent. Such a variant was found in our previous work 
p7j| , where a step factor: 

_ (prf 
Jtr — 




1 + {prf 

has been used for the truncation. The step-like truncation leads however to technical difficulties when 
used in a three-body calculation. Therefore, we present here some alternative variants of the triplet-even 
potential which still have an exponential truncation factor (1 — exp(— ar)) n with n > 3, for which the relative 
magnitude of the D-wave loop varies from 0.5 down to 0.02. It is important to note that the stronger the 
truncation of the tensor potential, i.e. the weaker the tensor interaction at small distances, the larger the 
3N binding energy (see the subsequent section). 

The dependence of all observables on the truncation factor of the central part V° PE is very weak. For 
example, a stronger truncation in Vq PE is easily compensated by a corresponding strengthening of the 
central attractive (Gaussian) part of the interaction. Therefore, in the present three-body calculations, we 
have used only one of the new variants for the singlet potential presented in Table 3 (namely, variant c). 
In addition, we have found that a fit to the experimental phase shifts in the interval 0-400 MeV does not 
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uniquely determine the width r\ and depth Vq of the deep potential in the central part of the NN interaction 
(35). Any variation of the width 77 can always be compensated by an appropriate variation of the depth Vo, 
without changing the quality of the phase shift fit. As a result, it is possible to choose the same width for 
the triplet 3 Si — 3 D\ (r/t) and the singlet 1 So — 1 D$ (r) 8 ) channels. Alternatively, it is possible to choose 
the same values of the depth Vq in triplet and singlet channels. In this case the width rj s and r\t would be 
different. 

We emphasize that the equality r\ s = rjt is impossible to reach in simple conventional potential models 

0. In other words, the width of the triplet potential is always different from that of the singlet potential. 
Both parameters (i.e. width and depth) of the conventional potential model are uniquely determined by the 
values of the low-energy effective range parameters, i.e. by the scattering length and effective range. It was 
recently found 55 that owing to some freedom in the choice of the three main parameters of the interaction, 

1. e. width 77, depth Vq, and the truncation radius of the OPE potential a, the values of the width parameters 
r\ and a 2 can be taken to be equal, i.e. r\ — a 2 . We then obtain the simplest two-parameter model of the 
Moscow potential, where only two parameters (width and depth) are fitted to the scattering length and 
effective range of the 1 Sq partial wave. However, unlike conventional force models (see ref. jjjj Chap. 2), we 
obtain with this simplest two-parameter model an excellent description of the 3 Si — 3 I?i and ^So-wave phase 
shifts in a large interval of 0-500 MeV instead of a good fit only in the range between 0-15 MeV typical 
for conventional models. Simultaneously, we obtain a good description of the main static properties of the 
deuteron |p5j| , including the quadrupole moment, charge r.m.s. radius and the asymptotic normalizations 
As and Ad- 

We attribute this success to the choice of the proper degrees of freedom for a potential description. Con- 
ventional models try to describe both, the short-range six-quark and the long-distance one-boson exchange 
aspects of the NN interaction. This leads to a complicated energy and momentum dependence of the re- 
sulting NN potential. We claim that the Moscow potential model is so simple, because we do not try to 
cast the six-quark aspects of the NN interaction into the same formalism as the asymptotic long-range part. 
The six-quark aspects in the NN sector are more adequately described by an orthogonality condition, rather 
than by a local repulsive NN potential. 

The parameters of the Moscow potential model are listed in Table 3. The deuteron wave function calculated 
with model B is shown in Fig. 4, and its static properties are listed in Table 4. In Fig. 5, the phase shifts 
of the even partial waves for variant (B) are compared with the experimental ones |57[ . The description of 
the NN phase shifts is for all variants listed in Table 3 quite reasonable. However, our description of the 
mixing parameter t\ especially for energies E > 400 MeV, is not yet satisfactory. The mixing parameter e\ 
is determined by the value of the matrix element s\Vt\^ d)- Our overestimation of t\ at higher energies 
is a consequence of the behavior of the 3 Si and 3 Z?i eigen functions at short distances. While the S-wave 
has a node and a rather pronounced loop, there is almost no loop in the D-w&vc. This leads to a rather 
sharp increase of e\ at higher energies. ^ In the more complete two-channel model, the short-range S-wave 
loop should be significantly reduced and we expect a flatter behavior of e± in better agreement with modern 
data. 



B. Exclusion of the deep- lying bound state 

The "extra" bound states in S- and P- waves have to be excluded if the potential is used in few-nucleon 
calculations. As was already discussed in Section 2, this is most conveniently done by means of the orthog- 
onalizing pseudo-potential (OPP): 

Vopp =XT, r = \<p f X<pf\ (40) 

with a large positive value of the strength constant A pj,^9|. In practice a value of A ~ 10 5 — 10 6 MeV is 
quite sufficient. In the case of coupled channels 3 Si — ^Di and 3 P2 — 3 F2 the function iff in (40) has two 
components. Therefore, it is necessary to use a two-channel matrix projector of the form: 



\ \<p 2 f >< <p)\ \tpj >< <p}\ j 



(41) 



where \(p^ > and \ip 2 > are the upper and lower entries of a two-component column vector, corresponding to 
the "extra" bound state. However, the matrix projector considerably complicates three-body calculations. 
Therefore, we use instead a one-component projector: 



J In contrast to this, in the first version of the Moscow potential model ^8|, there were two coinciding nodes and 
loops in both the S and D waves and a rather satisfactory fit to ei was obtained. 
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(42) 



acting on the "main" channel, i.e. on the channel possessing highest weight in the "extra" state. In other 
words, in 3 Si — 3 Di channels we use only an S-wave projector. In this case the function (pf is selected so 
that the inclusion of the above one-component projector into the Hamiltonian gives a result that is as close 
as possible to the action of the two- component eigenprojector IV Strictly speaking, the replacement of the 
matrix projector (41) by a single-channel projector (42) should cause some modification of the appropriate 
NN phase shifts. We have found, however, that this replacement influences the description of the phase 
shifts in coupled channels only very weakly. Moreover, employment of the simple Gaussian function with 
suitable values of r in (pf can ensure a quite reasonable and accurate description of scattering phase shifts. 

We will show, in particular, how the approximate S'-wave projector T = \(pf >< (p/\ acts in the coupled 
3 Si — 3 Z?i channels. The dependence of the first two energy levels on the orthogonalizing coupling constant A 
for the 3 Si — 3 Z?i channels is shown in Fig.6(a-b) for the potential (B). If the projector is absent, there are two 
bound states: the "extra" ground state with energy Eq = —599 MeV with a D-wave probability of 6.3%, and 
the second bound state describing the physical deuteron with E\ = —2.2245 MeV and a D-wave probability 
of Pd=5.75 %. If the exact two-component state ipQ is used in T\, the energy of the ground state Eq is 
shifted by an amount A when A is varied, whereas the deuteron energy E\ does not vary at all (by definition). 
The deviation of the curves Eq(X) and -Ei(A) from the straight line characterizes the difference between the 
application of the approximate (single-channel) S'-wave projector and the exact matrix eigenprojector IV 
As one can see, this difference is only important in the narrow range of values A = 685 — 690 MeV, where 
the energies Eq and E\ of both states come close to each other and the states are mixed strongly due to 
approximate orthogonality of (pf and ipi (their overlap being < <p>f\(pi >= 0.001). For A > 690 MeV, the 
deuteron state becomes the new ground state and its energy grows very slowly with A approaching the 
limiting value E(oo) = —2.2246 MeV. On the other hand, the original ground state energy Eo(X) continues 
to rise and finally transforms into a narrow resonance. In the case of the exact eigen projector, it would 
be transformed into a bound state embedded into the continuum with the positive energy Eq + A. This 
resonance strongly distorts the phase shifts only in a narrow energy range near E cm — Eq + A. But since the 
values for A used in the present three-nucleon calculation are of the order of A ~ 10 6 MeV, this distortion 
is far away from the physically relevant region. It is worth to note that if one takes only the S-component 
of the exact wave function ifQ as tpf in the single-channel projector (42), the result turns out to be much 

worse. In other words, the result of the application of the projector >< <Pq \, i.e. using the exact 

first component of two-component vector, differs much more from the matrix eigenprojector F2, than the 
approximate projector \ipf >< tpf\. This is explained by the fact that the orthogonality condition which is 
automatically satisfied for the eigenprojector fa is: 



From here it does not follow, in general, that < ip \ip\ >= 0. At the same time it is possible to make 
the overlap of the deuteron with some approximate function ipf arbitrarily small. Then the energy E\ will 

not depend on A, as in the case of the application of the exact eigenprojector Y^- However, it is impossible 
to choose a simple approximate function that will be orthogonal to the scattering wave functions of the 
continuous spectrum for all energies. Therefore, due to the employment of the approximate projector, the 
scattering phase shifts are inevitably distorted in some small energy interval. 

Bearing in mind the use of the potential in the three-nucleon system, where large values of the constant 
A are needed, we restrict ourselves to the approximate projector with a simple Gaussian form factor. Then 
the additional pseudo-potential V sep has the form of eq.(40) in all channels. For S- and P partial waves, 
the constants A have to be sufficiently large (~ 10 5 — 10 6 MeV) in order to ensure the exclusion of "extra" 
bound states, while their values in D-waves are determined by fitting the experimental phase shifts. The 
corresponding values of tq for the approximate projector are also given in Table 3. In the 3 Z?i channel and 
in all higher partial waves with I > 3, the term V sep is absent. 



Here, we discuss the predictions of the new versions of the Moscow model described in section 4. The 
3N calculations have been done by a variational method on a very large, nonminimal and nonorthogonal 
Gaussian basis [^2],[30| and are, in general, rather similar to the very accurate 3N-calculations of Kameyama 
et al. [ p3[ . With this basis the absolute accuracy for the eigenenergy is about 100 keV. The relative accuracy 
for the Coulomb displacement energy is even higher. 



< $ |$i >=< f i o 1) \fi ) > + < <^o 2 Vi 2) >= 0. 



(43) 



V. PROPERTIES OF THE THREE-NUCLEON BOUND STATE 
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As follows from many previous 3N bound state calculations, the total contribution of odd partial wave 
interactions to the 3N binding energy is small and does not exceed 0.2 MeV ||. We emphasize that odd 
partial waves have been included in our variational basis for the 3N bound state calculation but we did not 
take into account the contribution of the odd partial waves in the NN force. In the NN force, we include 
all even partial waves up to total orbital angular momentum L — 4. Our results are presented in Table 5, 
where we list several static properties of 3 H and 3 He such as binding energies E, root-mean-square charge 
radius r c h, the percentages of D- (Pd) and P- waves (Pp), and the Coulomb displacement energy Acw for 
the difference Eb( 3 H) — Es^He). In addition, we list the average kinetic energy T in the ground state, as 
well as certain characteristics of the 3 H charge form factor, i.e. the position of the first minimum q m im the 
second maximum q max and the value of \F c h\ in the second maximum F max . 

A. Dependence on the projection constants 

First, we discuss the dependence of the three- nucleon properties found with the above force model on 
the projection constants A. In Fig. 6(a-b) we show the dependence of the 3N ground state energy E 3 n on 
the value A( 3 Si) multiplying the projector V, projecting onto an "extra" bound state in the 3 Si channel. 
The two-body state energies E (X) and £a(A) shown in Fig. 6 were already discussed in Section 4. All 
other parameters of the potentials are fixed, including the constants A for the other channels. For A <680 
MeV the 3N-system remains unbound and the variational value E^n practically coincides with the energy 
of the two-body ground state E . When A is increased beyond 682 MeV, the system becomes bound with 
E-jn = — 6 MeV. If A is increased further the three-nucleon binding energy decreases very slowly approaching 
its limiting value £^(00) = —5.74 MeV, because these test calculations were done in a comparatively small 
basis. The saturation occurs for A = 10 5 — 10 6 MeV. For higher values of A the numerical results become 
unstable. The convergence of variational calculations with respect to A was investigated in some detail in 
Ref. The physical deuteron becomes the ground state of the two- nucleon subsystem for A >688 MeV 

when \E 3 n\ is already less than 6 MeV. 

These results clearly show that it is impossible to fit the binding energy of the three-body system to the 
experimental value by adjusting the parameter A. It is interesting to note that the system becomes bound 
when the repulsion in pair subsystems increases^. The reason of this interesting phenomenon is the following. 
When the projection constant in the singlet channel is chosen as e.g. A s = 10 6 MeV and the one in the triplet 
channels as e.g. At << 10 6 MeV (and varying), the almost complete space symmetry of the three-nucleon 
system is broken, and the weight of wave function components with mixed orbital permutational symmetry 
[21] x is increased. 

On the other hand, if the constants A( 3 Si) and A( 1 S'o) evolve simultaneously, we find a different behavior 
-E"3Jv(A) that is shown by the dashed line in Fig. 6b. Now, the 3N system remains bound for all A. However, 
for the allowed values of A that is when both "extra" bound states in the triplet and singlet channels are 
pushed sufficiently high above the physical state, the three-nucleon binding energy \Esn\ appears at a value 
less than 7 MeV. So, in this case it is also impossible to fit the three-nucleon binding energy by shifting 
the "extra" two-body state to the positive energy region. Therefore, all results for the three-nucleon system 
discussed below have been calculated with sufficiently large values X s and A*(~ 10 5 — 10 6 MeV) for which 
complete saturation has been reached, and the results are, in some sense, A-independent. 

B. Binding energy of 3 H and 3 He. 

All variants of the model studied here differ only in the form of the truncation factor for the OPEP-tensor 
force at small distances and lead to overall similar results (see tables 3,4,5) except for one previous version 
of our model |5(J (variant A in Tables 3-5). In the course of the calculations we have found that the stronger 
cut-off in OPEP tensor force, the larger the 3N binding energy (see the Tables 3 and 5). In this respect our 
results are very similar to those found by Sasakawa and Ishikawa pij . In their 3N studies with the Reid 
Soft-Core (RSC) NN potential using various truncation factors of the tensor force they found that when the 
tensor force is cut off more strongly at short distances (and the central attractive part is correspondingly 
increased in order to keep the deuteron binding energy invariant), the 3N binding rises as well. 

However, there is an important difference between the present results and those in Ref. In our 

case, the partial replacement of the short-range tensor force by a central force does not lead to a noticeable 



Certainly, the absolute value of three-body binding energy decreases in the process, however the energy of the 
two-body threshold decreases still more quickly. 
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reduction of the tensor mixing parameter e\. In fact as the truncation is increased, the mixing parameter 
does not tend to zero, but, instead of this, it rises more sharply for larger n. This is related to the essentially 
different character of interference between tensor and central force (both in 2N and 3N sy stems) in our case 



as compared to the traditional nuclear force models (see the last paragraph in Sec . [IV A| ). However, we do 
not expect drastic effects of our t\ description on the 3N binding energy, at least not within the corridor 
given by the different versions of our model. In fact, the difference in ei values between the RSC and the 
Paris potential model is of the same size as that between the present version and experiment j6^ ]. 

The 3 H binding energies obtained in the present study (~6.05 MeV) for the two most realistic variants of 
the potential (B and C) are much higher than the value of ^4.5 MeV obtained by Hahn et al. Q] in their 
first 3N calculation with the earlier version of the Moscow potential. This large difference is caused by the 
following reasons: 

(i) In the Hahn et al. |Q calculation and also in the early 1983 version of the Moscow potential an old 
version of the singlet 1 So potential fitted to the 1983 phase shift analysis [[H was employed. Our new singlet 
potential is fitted to the SAID database of 1995 |57[] and resulted in some increased 3N binding (about ~0.2 
MeV). 

(ii) Another, more important reason for the present improvement of the 3N results is the use of a stronger 
truncation of the short-range OPEP tensor force and the respective reduction of Pd in the deuteron and 
3N systems - see the versions B, C and D in Table 4. The sharper truncation removes or supresses to a 
significant degree the inner node and the respective loop in the D-wetve at short distances. 

(iii) A very restricted configuration space of only 5 basic channels 3 Si — 3 D\ and 1 Sq was employed in the 
Faddeev calculation of Hahn et al. [Q and the contribution of all I >2 higher partial waves was neglected. 
The latter contribution are quite important in our case especially in the Faddeev like 3N-approach. 

We have previously emphasized the extremely important contribution of higher partial waves as manifested 
by a sharp increase of the average 3N kinetic energy. Compare the third column in Table 5 with the 
corresponding results |L0]]64|] for traditional force models, in particular for model B: 

(Paris) = 42.6 MeV; E 3 k g(RSC) = 49.3 MeV; 



while Elg (Moscow) ~ 150 MeV, 

is more than in three times higher. It is worth to recall here that the modest 20% increase in the average 
kinetic energy in the RSC potential as compared with the Paris potential, results already in an enhancement 
of the contribution of higher partial waves to the binding energy by as much as a factor of 2.5. Thus, one can 
easily imagine how large a contribution from higher partial waves may result in our case. A more detailed 
analysis of these results will be carried out elsewhere. 

We must keep in mind that, contrary to conventional force models, our NN-potential corresponds only to 
a one-channel approximation of the complete two-phase model, that is by definition of the NN potential, 
only the effective NN channel with unexcited nucleons. As was argued above, the effective channel potential 
is phase-shift equivalent to the initial two-phase interaction model of hybrid type (see eqs. (19) and (24)) 
and the replacement of the local repulsive core by a two-channel interaction model (e.g. as in case of QCB 
model by Simonov |ll|) will lead to some additional binding of the order of some ~ 0.8 ~ 0.9 MeV |6j|, 
simply because of the presence of the second channel. There is an additional reason for extra-binding in 3N- 
and generally in many-nucleon systems if the full two-phase interaction model is considered. One has to take 
into account the interaction between the six-quark bag and the spectator nucleon. This dynamical 6q-N 
channel coupling can be considered as some kind of three-body force that will also increase the 3N-binding 
energy. We expect, that these effects will mainly change the absolute values of the energies while the relative 
energies, which will be discussed in the following, will be influenced only weakly. This is because the weights 
of the 6g-bag components in 3 H and 3 He are small and have almost equal magnititudc. 



C. Coulomb displacement energy, charge radii and form factors. 

As follows from the results displayed in Table 5, the Coulomb displacement energy: 

Acoui = E B { 3 H)-E B ( 3 He) 
(for variants B and C with n—5 and 7 for the OPEP tensor force truncation) equals: 

Acoui 670 keV 

for E B ( 3 H)=6.05 MeV. For a more careful test of the model, it is important, however, to find the value of 
Acoui when E B (^H) is scaled to its experimental value —8.48 MeV since the larger the 3N binding energy 
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the smaller the r.m.s. charge radius and the higher the Coulomb energy of 3 He. Thus, we slightly increase 
(by ~ 2%) the strength of the central (triplet and singlet) potential well to fit the experimental value of 
Eb( 3 H), and for this case we analyze the Coulomb energy of 3 He and the r.m.s. charge radii of 3 H and 3 Hc 
with our model (see the two last rows in Table 5). 

The extrapolated Coulomb displacement energy turns out to be: 

A^iT = 737 keV 

which is only by 27 keV less than the experimental value. Here, we note in passing that in the full two-phase 
model calculation the magnitude of Acoui would increase further because the second channel comprises the 
more tightly packed 6q + N configuration. 

The most likely reason for this improvement in the prediction of A^f^f is the appearance of the inner 
radial nodes and corresponding short-range loops along every interparticle coordinate (and also along 
both Jacobi coordinates) in our 3N wave functions. The inner loops along the r»j coordinates have only a 
minor influence on the r.m.s. radii of the charge or matter distribution because the latter are determined 
mainly by one-particle coordinates (as measured from the center of mass). Contrary to this, the Coulomb 
potential between two protons in 3 He is maximal just at zero interparticle separation where we have some 
enhancement of particle density due to the presence of the inner loops in our wave functions. 

This naturally explains why conventional force models leading to a short-range suppression of the wave 
functions along every interparticle distance miss some 100 keV in the Coulomb energy of 3 He, whereas our 
model, due to the inner loops leads to a prediction that is close to the experimental value. Although our 
result is still qualitative in the sense that it involves an extrapolation to the empirical binding energy, the 
observation that the inner wave function loops arc extremely important for the Coulomb displacement energy 
is likely to survive. Moreover, if the force model presented here is corroborated in subsequent studies, it 
would suggest that inner loops are ubiquitous in nuclear wave functions. 

The recent quantitative explanation of the Coulomb displacement energy in the literature |j66f using a force 
model with a repulsive core is based on the assumption of a rather strong charge symmetry breaking (CSB) 
term Vcsb (difference between np and nn strong interaction). With Vcsb added to the the Argonne Vis 
potential, one can almost completely explain the 100 keV difference in A cou i between 3 He and 3 H. However, 
the same Vcsb leads to a strong overestimation of the Coulomb energy difference in the 6 Be — 6 He system. 
Therefore, it is conceivable that the contribution of the charge symmetry breaking term Vcs b in Ref . |6(| is 
actually smaller. In any case, it is evident from the discussion above that the assumption of a strong charge 
symmmetry breaking interaction Vcsb is not the only mechanism that can explain the 100 keV Coulomb 
energy difference in the 3N system. 

The values of the extrapolated charge radii of 3 H and 3 He (see last row in Table 5) turn out to be 
larger than the corresponding experimental values. However, it should be emphasized once again that the 
3N-calculations with one-component force model presented here, i.e. carried out only in the NN (or 3N-) 
sector, corresponds just to the projection of the total two-phase model wave function onto the pure 3N- 
component. The two-phase solution includes also compact 6q + N components with a sizeable probability, 
which may considerably reduce the 3 A charge radius. Similarly, the full two- phase model will lead to 
important modifications in the charge form factors of 3 H and 3 He. Our results for the 3 H charge form 
factors calculated with the one-phase Moscow potential are displayed in Fig. 7 together with the predictions 
for the extrapolated (to experimental Eb( 3 H)) solution. We see that the behavior of the form factor for 
q 2 > 14 fm~ 2 is rather similar to predictions calculated with a conventional force model such as the RSC 
potential. However, we expect appreciable modifications for q 2 > 14 fm -2 for the original two-phase model, 



/p-pr 

in whic h the s ign of the 6q + N component is opposite to that of the main 3 A-component |13| . According 



to refs. 1 65 67]] this fact should considerably improve the description of the 3 A charge form factors in the 
secondary maximum. 

Summarizing we conclude that the description of both the 3A binding energy and some important 3A 
observables will presumably improve when passing from a one-phase potential model to a two-phase model 
due to the explicit appearance of the 6q+N component in the three-body solution. 



VI. CONCLUSION 



The force model presented in this paper -referred to as M-type model- differs in a few important aspects 
from traditional NN interaction models currently in use. 

First, M-type models necessarily incorporate additional orthogonality condition(s) with respect to cer- 
tain nodeless functions ipQ. The orthogonality condition can be interpreted as the projection of six-quark 
compound states, ipa with maximal space (or permutational) symmetry onto the asymptotic NN-channel. 
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Second, M-type models are characterized by the presence of a deep attractive well while a part of the 
short-range repulsive core (or the whole core) is replaced by appropriate orthogonality conditions^} The 
consequences of the different off-shell behavior of conventional and M-type NN force models can be tested 
in NN bremsstrahlung 

Finally, the models presented here have a physical interpretation that is rather different from that of 
conventional force models. A standard force model aims at describing the NN interaction in the total 
configuration space in the form of some potential. Contrary to this, we start from the assumption that 
certain six-quark compound states of maximal permutational symmetry are present in the lowest partial 
waves. The latter have only very small projections onto the non-antisymmetrized asymptotic NN channel. 
However, at the same time they are strongly coupled by the six-quark Hamiltonian with the clusterized 
NN channel at small distances r <1 fm. As a result, these compound states cannot be described by any 
reasonable NN potential. These states can only be described by a very complicated nonlocal and energy- 
dependent NN interaction operator. However, even when utilizing such complicated interaction operators 
we are forced to introduce many adjustable or free parameters to fit the data. Thus, due to the formation 
of the above compound six-quark states, the whole system cannot be described by any simple NN potential. 
In order to get an efficient potential description, it is necessary to remove the above nonpotential pieces 
from the full interaction operator. This can be done with the help of the two-phase orthogonalized model 
described in Sect. 2. The exclusion of six-quark compound states can be conveniently accomplished by the 
well known Feshbach formalism. This eventually leads to the Moscow type models with its two orthogonal 
channels. 

Our approach makes it possible to describe rather accurately both NN phase shifts up to 500 MeV, and the 
deuteron structure with a truncated one-pion exchange potential together with a simple deep static potential. 
Altogether we use 6 parameters with clear physical meaning. This evidently confirms the usefulness of the 
present two-phase approach. Other hybrid models which combine both quark- and meson exchange degrees 
of freedom make use of nonorthogonal quark- and meson exchange channels, i.e. some mixture of both. As 
a result, they contain nontrivial energy dependence and nonlocalities. In addition, most hybrid models do 
not offer a microscopic interpretation of the NN channel. 

A clear separation of the nonpotential pieces of the NN interaction and the subsequent parametrization of 
the rest in an orthogonal subspace is the main physical idea underlying our approach. Therefore, the success 
of the model may be taken as some evidence for the formation of six-quark compound states (i.e. dibaryons) 
at short-range in low partial waves. As a result of the elimination of 6q-compound states, the 3N-, 4N- and 
other many-nucleon systems are underbound. 

In this work we have used only the one-phase effective two-body potential component of the complete 
two-phase model. In the future, it would be interesting to take the next step and to incorporate the second, 
i.e. 6g+N component in few- and many-body calculations. In addition, it would be interesting to include 
vector mesons into the theoretical description. Due to the additional orthogonality condition, their effect 
will not be as large as in standard force models, where the wNN coupling constant is in contradiction with 
the SU(3) prediction H. In our case we are able to keep the wNN coupling constant as low as SU(3) predicts. 



Acknowledgments. 

We are thankful to many of our colleagues in Moscow and Tubingen for long term fruitful discussions 
on the topics of the present study, especially to Profs. E. W. Schmid, V.G.Neudatchin, and Dr. I. T. 
Obukhovsky. One of the authors (V.I.K.) is very grateful to Profs. Steven Moszkowski and Gerry Brown 
for many useful comments on our approach. He is also grateful to Prof. C.N.Yang, Head of Theoretical 
Physics Institute at Stony Brook, for his invitation to the Institute. The Russian authors thank the Russian 
Foundation for Fundamental Research (grant No. 96-02-18071) and the Deutsche Forschungsgemeinschaft 
(grant No. Fa-67/20-1) for partial financial support. 



7 The M-type models presented here are differ from the OCM approach by Saito |Q in a few essential aspects: 
our model is based on another, more convenient mathematical formalism that is used to incorporate the additional 
orthogonality conditions. Second, the physical meaning of the "forbidden" subspace is quite different in our case 
because we start from a two-phase interaction model. Third, we generally combine the (reduced) repulsive core and 
the additional orthogonality condition. 
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Figure captions 



Fig. 1. The scattering wave functions at Ei a b—1 MeV for three S'-wave phase shift equivalent potentials: 
initial local MT-V (curve 1), the intermediate nonlocal potential with ro=0.2 fm (curve 2), and the purely 
attractive potential with the projector close to the eigen projector with ro=0.4 fm (curve 3). 

Fig. 2. NN S-wave phase shifts for center of mass energies between 0-400 MeV for the initial MT-V 
potential (dashed line) and for two modified potentials V'(r) fitted to S- and .D-wave phase shifts (see Table 
2): (i) with ro=0.15 fm (solid line) and (ii) with ro=0.2475 fm (dotted line). 

Fig. 3. D-wave phase shifts for the same three variants of the interaction potential as in Fig. 2. 

Fig. 4. The deuteron S-wave (1=0) and D-wave (1=2) functions for model B. The deuteron wave functions 
calculated with the Paris potential B are shown for comparison. 

Fig. 5. The even parity phase shifts for the new version of the NN Moscow potential (variant B) are 
compared with the data of the energy-dependent phase-shift analysis by Arndt et al. |57j] (triangles, dotted 
line for £i). For E\ we also show the Nijmegen phase-shift analysis PWA93 (5£| (dashed line) and the single 
energy phase-shift analysis by Arndt et al. |57j] (circles). 

Fig. 6. (a) Effect of the approximate S'-wave projector on the discrete spectra of 2N and 3N systems. Eo 
and Ei are the energies of the first two 2N-states, E^n is the energy of 3N ground state, A is the 3 Si-wave 
projector coupling constant, (b) The dotted line shows the dependence of the 3N-bound state on A for the 
case A s = At. 

Fig. 7. The charge form factor of 3 H for different versions of the Moscow NN potential: the previous 
version (A) (short-dashed line), the present version (variant B) (dashed line), and the potential extrapolated 
to E exp ( 3 H) (solid line). Also are shown are the previous (□) and the new |6^] (A) experimental data. 
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TABLE I. The parameters of the modified potentials reproducing the S-wave phase shifts for the initial MT-V 
potential with different values of the "projector radius" ro and properties of the wave function. 



ro, [fm] 


g' R [MeV] 


9a [MeV] 


r-node [fm] 




(1 MeV) 


per point 


e [MeV] 


eo[MeV] 


0.00 


1458.1 


-578.1 


0.00 






-0.4136 




0.15 


992.0 


-511.0 


0.25 


0.03 


0.40 


-2.51 


-0.385 


0.18 


673.8 


-468.6 


0.28 


0.06 


1.07 


-7.55 


-0.370 


0.20 


433.0 


-438.0 


0.30 


0.07 


1.70 


-36.06 


-0.364 


0.25 


-5.0 


-389.0 


0.35 


0.16 


2.50 


-430.30 


-0.399 


0.30 


-63.6 


-403.0 


0.40 


0.28 


1.60 


-683.70 


-0.611 


0.40 


-106.2 


-494.0 


0.50 


0.47 


0.05 


-528.00 


-0.919 


0.50 


-287.2 


-582.5 


0.60 


0.63 


0.10 


-378.80 


-0.880 



TABLE II. Variants of the MT-V potential with different values of the "projector radius" ro, fitted to S- and 
D-wave phase shifts of the original MT-V potential and the two-body bound state energy £o=-0.4136 MeV, together 
with the corresponding results for the three-body calculations. 



ro 


9'r 


9'a 


E 3 


T 3 


%S 


%D 


Si rprep 


Qmin 


Qmax 


+ max 


fm 


MeV 


MeV 


MeV 


MeV 






MeV 


fm" 2 


fm" 2 




MT-V 


1458.05 


-578.09 


-8.25 


30.6 


99.0 


0.75 


-7.52 


19 


24.0 


4.8-10" 4 


0.15 


999.36 


-514.00 


-8.42 


40.3 


98.5 


1.1 


-6.73 


19 


24.0 


5.4-10" 4 


0.18 


800.00 


-497.76 


-8.45 


44.3 


98.0 


1.3 


-6.11 


18 


23.5 


6.0-10 -4 


0.20 


445.80 


-443.30 


-8.82 


51.8 


97.8 


1.7 


-5.96 


19 


25.0 


6.0-10" 4 


0.25 


0.00 


-388.67 


-8.99 


74.6 


97.0 


2.0 


-5.36 


20 


25.0 


5.9-10" 4 


0.30 


-71.50 


-394.50 


-8.02 


115.0 


95.0 


3.0 


-4.20 


21 


27.0 


4.0 ■ 10~ 4 


0.40 


-69.87 


-464.90 


-7.02 


189.0 


89.0 


5.8 


-3.20 


24 


31.0 


1.7-10 -4 



TABLE III. Parameters of the Moscow NN potentials. 



variant 




parameter values 








parameter values of ro [fm] 


A [MeV] 




n 


Vo [MeV] 


V [fm-'] 


a[fnr L ] 










for particular channels 




triplet channels 










D 2 




A 


3 


-466.7 


1.600 


3.000 


0.448 




oo 


0.4600 


251.0 


0.300 


oo 


B 


3 


-1329.0 


2.296 


1.884 


0.445 




oo 


0.4061 


427.2 


0.276 


oo 


C 


7 


-1459.2 


2.372 


2.610 


0.454 




oo 


0.4199 


353.8 


0.289 


oo 


singlet channels 




L S 




L D 2 




a 




-1106.2 


1.600 


3.000 


0.4998 




oo 


0.4472 


229.0 






b 




-1220.0 


1.753 


1.884 


0.4815 




oo 


0.4103 


303.7 






c 




-1222.0 


1.738 


1.000 


0.4825 




oo 


0.4071 


321.8 







TABLE IV. Deuteron predictions for the different versions of the triplet-even Moscow NN potentials. Di oop is 
the amplitude of the D-wave loop. The deuteron matter radius is defined as r 2 ^ = 1/4 J" 00 drr 2 (u 2 (r) + w 2 (r)). For 
the most recent value of the deuteron matter radius see Ref. [Q. Here, r\ is the asymptotic D/S ratio. In impulse 
approximation (model B) the deuteron magnetic moment is fid = 0.847/ijv, the quadrupole moment is Qd = 0.271 
fm 2 , and the charge radius is r c h = 2.112 fm. 



Variant 


x 2 


Dloop 


E d [MeV] 


Pd% 


n 


r m [fm] 


A 


118 


0.50 


-2.2244 


6.59 


0.0267 


1.96 


B 


56 


0.30 


-2.2246 


5.75 


0.0258 


1.95 


C 


126 


0.02 


-2.2246 


6.14 


0.0262 


1.94 



22 



TABLE V. Three-nucleon properties with the A, B, C-versions of the Moscow potential. 



variant of 


E( 6 H) 


T( A H) 


Pd 


Pp 


r ch ( A H) 


A.Coul 


rchCHe) 


•2 


'2 




potential 


MeV 


MeV 


% 


% 


fm 


MeV 


fm 


far 2 


fm" 2 




A 


-5.07 


143.7 


7.8 


0.03 


2.15 


0.612 


2.41 


13 


18 


1.2-10~ ;i 


B 


-6.03 


162.4 


7.4 


0.03 


2.01 


0.670 


2.24 


16 


22 


9.6-10~ 4 


C 


-6.03 


162.8 


8.1 


0.05 


2.01 


0.667 


2.25 


16 


22 


9.6-10" 4 


extrapol. 


-8.48 


175.1 


8.5 


0.06 


1.88 


0.737 


2.07 


17 


23 


1.0 • 10" 13 


to E exp 






















experiment 


-8.48 








1.76 


0.764 


1.96 


13 


18 


4 • \Q- A 
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. I — I 

o 

c 
3 



CD 
> 



1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 
-0.25 
-0.50 
-0.75 
-1.00 
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